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fvq ' Abstract. A recent theory is reviewed for the saniple-to-sample fluctuations 

in the critical current of a Josephson junction consisting of a disordered point 
contact or microbridge. The theory is based on a relation between the super- 
'Trt ' current and the scattering matrix in the normal state. The root-mean-square 

(~| ' amplitude rms/c of the critical current /c at zero temperature is given by 

rms/c ~ eAo/h, up to a numerical coefficient of order unity (Aq is the en- 

l5 ' ergy gap). This is the superconducting analogue of "Universal Conductance 

Fluctuations" in the normal state. The theory can also be applied to a ballis- 
tic point contact, where it yields the analogue of the quantized conductance, 

C^ , and to a quantum dot, where it describes supercurrent resonances. All three 

phenomena provide a measurement of the supercurrent unit cAq/U, and are 
"universal" through the absence of a dependence on junction parameters. 

a, 

2 ! 1 Introduction 

Nanostructures combining semiconducting and superconducting elements form 
\N ■ a new class of systems in which to search for mesoscopic phenomena. The corre- 

J^ I spondence between transport of normal electrons and transport of Bogoliubov 

pg ■ quasiparticles (being the elementary excitations of the superconductor) serves 

as a useful guide in the search. The appearance of a new length scale — the 
^^D ' superconducting coherence length i^ — complicates the correspondence in an 

^P I interesting way, by introducing two qualitatively different regimes: The short- 

junction regime (junction length L <C ^) and the long-junction regime {L ^ £^). 
These two regimes appear over and above the transport regimes which depend 
j^ ' on the relative magnitude of junction length and mean free path I: The ballistic 

regime {I 3> L) and the dijjusive regime {I ^ L). A third transport regime in 
the normal state, that of resonant tunneling, is also subdivided into two new 
' O I regimes, distinguished by the relative magnitude of ^ and the characteristic 

^ ■ length vpTj-cs, with Ties the lifetime of the resonant state in the junction and vp 

fj I the Fermi velocity. 

Our own interest in this field has focused on a set of three phenomena which 
provide, through the Josephson effect, a measurement of the supercurrent unit 



eAo/?i (Ao being the superconducting energy gap). This unit of current plays 
the role of the conductance quantum e^ /h in the normal state. The three 
phenomena are: 

1. Discretization of the critical current of a ballistic point contact [1]. 

2. Resonant Josephson current through a quantum dot [2, 3]. 

3. Mesoscopic supercurrent fluctuations in a diffusive point contact [4]. 

These effects belong, respectively, to the ballistic, resonant-tunneling, and dif- 
fusive transport regimes, and within each regime to the short-junction limit. 

Effect number f is the analogue of the quantized conductance of a quantum 
point contact [5]-[7]: The critical current /c of a short and narrow constriction 
in a superconductor increases stepwise as a function of the constriction width, 
with step height eAo/?i independent of the properties of the junction [1], 

The integer N is the number of transverse modes at the Fermi level which can 
propagate through the constriction. The short-junction limit i <C ^o (where 
^0 = hvp/nAo is the ballistic coherence length) is essential: Furusaki et al. [8] 
studied quantum size effects on the critical current in the long-junction limit 
and found a geometry-dependent behavior instead of Eq. (f ). 
Effect number 2 is described by the formula [3] 

^ e ApFp 

hAo + To ^ ' 

for the critical current through a bound state at the Fermi level which is coupled 
with equal tunnel rate To/h to two bulk superconductors. Since t^^s ~ ^/'^Tq, 
the "short-junction" criterion ^g 3> wpTres is equivalent to Fq 3> Aq. If Fq 3> 
Ao, the critical current on resonance becomes [2, 3] 

ir = ^, (3) 

which no longer depends on the tunnel rate. 

Effect number 3 is the analogue of "Universal Conductance Fluctuations" 
in disordered normal metals [9]-[ll]. The sample-to-sample fluctuations in the 
critical current of a disordered point contact have root-mean-square value [4] 

rms/c~— -, (4) 

up to a numerical coefficient of order unity. These mesoscopic fluctuations are 
universal in the sense that they do not depend on the size of the junction or 



Mesoscopic Joscphson Effects 



ballistic 
res. tunneling 
diffusive 



short-junction long-junction 



/c = NeAo/h geometry dep. 

/™^ = eAo/h 4^°^ = eVo/h 

rms/c ~ eAo/?i rms/c — eD/L^ 



Er. 



hvp/L 

To 
hD/L^ 



Table 1: Summary of the mesoscopic Josephson effects described in the text, 
in the ballistic, resonant-tunneling, and diffusive transport regimes. Universal 
(junction-independent) behavior is obtained in the short-junction limit, which 
in terms of the correlation energy Ec is given by Aq <C Ec ■ 



on the degree of disorder, as long as the criteria I <ti L <ti Nl and L <C C 
for the diffusive, short-junction regime are satisfied. (Here ^ = (CoO^^^ is the 
diffusive coherence length.) The short-junction limit L <C ^ is again essential 
for universality. The opposite long-junction limit was considered previously 
by Al'tshuler and Spivak [12], who calculated the supercurrent fluctuations of 
an SNS junction in the limit L ^ ^ (S=superconductor, N=normal metal, 
L=separation of the SN interfaces). Their result rms/c — ev^l/L'^ depends on 
both I and L and is therefore not universal in the sense of UCF. 

We have summarized the three universal short-junction Josephson effects 
in Table 1, together with their non- universal long-junction counterparts. We 
have introduced the diffusion constant D c^ v^l and the correlation energy (or 
Thouless energy) Ec — h/rdwcW, where Tdweii is the dwell time in the junction. 
The distinction "short-junction" versus "long-junction" may alternatively be 
given in terms of the relative magnitude of Ec and Aq: In the short-junction 
regime Aq ^ Ec, while in the long-junction regime Aq 3> -Be- The universality 
is lost in the long-junction regime because Ec (which depends on junction 
parameters) becomes the characteristic energy for the Josephson effect instead 
of Ag. All these results hold in the zero-temperature limit, or more precisely 
forfcBr<min(Ao,£:c). 

In the present paper we review a scattering theory for the Josephson effect, 
which provides a unified treatment of the three universal (short-junction) phe- 
nomena. In Sec. 2 we show, following Ref. [13], how the supercurrent can be 
obtained directly from the quasiparticle excitation spectrum of the Josephson 
junction. This approach is equivalent to the usual method which starts from 
the finite-temperature Greens function, but is more convenient in the short- 
junction regime, where — as we shall see — the excitation spectrum has a 
particularly simple form. In Sec. 3 we follow Ref. [4] in relating the excita- 
tion spectrum to the normal-state scattering matrix of the junction. In Sec. 
4 we take the short-junction limit and obtain a simple relation between the 
supercurrent and the eigenvalues of the transmission matrix product tt^ . This 
relation plays the role for the Joscphson effect of the Landauer formula for 



the conductance [7]. We emphasize that the transmission matrix t refers to 
normal electrons. Other scattering approaches to the Josephson effect [f4, f 5] 
relate the supercurrent to the scattering matrix for Bogoliubov quasiparticles. 
The usefulness of the present method is that one can use directly the many 
properties of t known from normal-state transport theory. Indeed, as we shall 
see in Sec. 5, the three universal Josephson effects follow as special cases of the 
general scattering formula of Sec. 4. 

2 Supercurrent from Excitation Spectrum 

The quasiparticle excitation spectrum of the Josephson junction consists of the 
positive eigenvalues of the Bogoliubov-de Gennes (BdG) equation [16]. The 
BdG equation has the form of two Schrodinger equations for electron and hole 
wavefunctions u(r) and v(r), coupled by the pair potential A(r): 

?4)(:)-(:)^ 

Here Hq — p'^ / 2m+V (r) — Ep is the single-electron Hamiltonian, containing an 
electrostatic potential V{r) (zero vector potential is assumed). The excitation 
energy e is measured relative to the Fermi energy Ep. 

In a uniform system with A(r) = Age"^, V{r) = 0, the solution of the BdG 
equation is 



e 



{{h^ky2m-EFf + Alf\ 
u(r) - V-i/2(2e)-i/V^/2(£ + ?i2fcV2m-i?F)'^'c'''", 
v(r) = V-i/2(2e)-i/V^/2(£-?i2fcV2m + i?F)'^%'''-'". (6) 

The eigenfunction is normalized to unit probability in a volume V, 

dr {\uf + \vf) ^ 1. (7) 



The excitation spectrum is continuous, with excitation gap Aq. The eigenfunc- 
tions (u, v) are plane waves characterized by a wavevector k. The coefficients 
of the plane waves are the two coherence factors of the BCS theory. This sim- 
ple excitation spectrum is modified by the presence of the Josephson junction. 
The spectrum acquires a discrete part due to the non- uniformities in A(r) near 
the junction. The discrete spectrum corresponds to bound states in the gap 
(0 < e < Ao), localized within a coherence length from the junction. In ad- 
dition, the continuous spectrum is modified. As we shall show in Sec. 3, in 
the short-junction limit the supercurrent is entirely determined by the discrete 
part of the excitation spectrum. 



Consider an SNS junction with normal region at \x\ < L/2. Let A(r) -^ 
Aoe^"^/^ for x -^ ±00. To determine A(r) near the junction one has to solve 
the self-consistency equation [16] 

A(r) = |5(r)|^v*(r)u(r)[l-2/(e)], (8) 

£>0 

where the sum is over all states with positive eigenvalue, and /(e) = [1 + 
exp(£/fcBT)]~"'^ is the Fermi function. The coefRcient g is the interaction con- 
stant of the BCS theory of superconductivity. At an SN interface, g drops 
abruptly (over atomic distances) to zero. (We assume non-interacting electrons 
in the normal region.) Therefore, A(r) = for |a;| < L/2 regardless of Eq. (8). 
At the superconducting side of the SN interface, A(r) recovers its bulk value 
Ap only at some distance from the interface (proximity effect). There exists a 
class of Josephson junctions where the suppression of A(r) on approaching the 
SN interface can be neglected, and one can use the step-function model 



Aoe'<^/2 


if a; < -L/2. 





if \x\ < L/2, 


Aoe-'"^/2 


iix> L/2. 



A(r) = <^ if |a;| < L/2, (9) 



The step-function pair potential is also referred to in the literature as a "rigid 
boundary-condition". Likharev [17] discusses in detail the conditions under 
which this model is valid: 

1. If the width W of the junction is small compared to the coherence length, 
the non- uniformities in A(r) extend only over a distance of order W from 
the junction (because of "geometrical dilution" of the influence of the 
narrow junction in the wide superconductor). Since non-uniformities on 
length scales <C C do not affect the dynamics of the quasiparticles, these 
can be neglected and the step- function model holds. Constrictions in the 
short-junction limit belong in general to this class of junctions. Note 
that, if the length and width of the junction arc <C £,, it is irrelevant 
whether it is made from a normal metal or from a superconductor. In 
the literature, such a junction is referred to as an ScS junction, without 
specifying whether the constriction (c) is N or S. 

2. Alternatively, the step-function model holds if the resistivity of the junc- 
tion region is much bigger than the resistivity of the bulk superconductor. 
This condition is formulated more precisely by Kupriyanov et al. [18], 
who have also studied deviations from the ideal step-function model due 
to the proximity effect and due to current-induced suppression of the 
pair potential. A superconductor — semiconductor — superconductor 
junction is typically in this second category. 



In equilibrium, a phase difference (f> in the pair potential induces a sta- 
tionary current / through the junction (DC Josephson effect). The current- 
phase relationship /(0) is 27r-periodic and antisymmetric [7(0 + 27r) = 1(4'), 
I{—4>) = —/((/>)]. The maximum I^ = max 7(0) is known as the critical current. 
There is a thermodynamic relation 

. '^e dF 

between the equilibrium current and the derivative of the free energy F with 
respect to the phase difference. To apply this relation we need to know how 
to obtain F from the BdG equation. The required formula was derived by 
Bardeen et al. [19] from the Greens function expression for F. An alternative 
derivation, directly from the BdG equation, has been given in Ref. [13]. The 
result is 

F = -2kBTj2^n[2cosh{e/2kBT)]+ fdr\A\^/\g\+Trno. (H) 

e>0 '' 

The first term in Eq. (11) (the sum over e) can be formally interpreted as the 
free energy of non- interacting electrons, all of one single spin, in a "semiconduc- 
tor" with Fermi level halfway between the "conduction band" (positive e) and 
the "valence band" (negative e) . This semiconductor model of a superconduc- 
tor appeals to intuition, but does not give the free energy correctly. The second 
term in Eq. (11) corrects for a double-counting of the interaction energy in the 
semiconductor model. The third term TrTio (i-c the sum of the single-electron 
eigenenergies) cancels a divergence at large e of the series in the first term. 
From Eqs. (10) and (11) one obtains [13] 

/ = -|^tanh(e,/2fcBr)^ 
p 



2e, 



2kBT rdeln[2cosh(£/2fcBT)]^ + ^-^ /dr |A|Vl5Ul2) 
Jao "<? Ti d<j) J 



where we have rewritten X]e>o ^^ ^ sum over the discrete positive eigenvalues 
Ep (p = 1,2,...), and an integration over the continuous spectrum with density 
of states p{e). The term TrTYg in Eq. (11) does not depend on 0, and therefore 
does not contribute to /. The spatial integral of |A|^/|(7| does contribute in 
general. If the step-function model for A(r) holds, however, |A| is independent 
oicf) so that this contribution can be disregarded. A calculation of the Josephson 
current from Eq. (12) then requires only knowledge of the eigenvalues, not of 
the eigenfunctions. 




Figure 1: Superconductor — normal- metal — superconductor Josephson junc- 
tion containing a disordered normal region. Scattering states in the normal and 
superconducting leads are indicated schematically. Normal scattering (disor- 
dered region) and Andreev reflection (SN interface) are spatially separated. 



3 Excitation spectrum from scattering matrix 

Following Rcf. [4] , we now show how to relate the excitation spectrum of Bo- 
goliubov quasiparticles to the scattering matrix of normal electrons. The su- 
percurrcnt then follows from Eq. (12). The model considered is illustrated in 
Fig. 1. It consists of a disordered normal region between two superconducting 
regions Si and S2. The disordered region may or may not contain a geomet- 
rical constriction. To obtain a well-defined scattering problem we insert ideal 
(impurity-free) normal leads Ni and N2 to the left and right of the disordered 
region. The leads should be long compared to the Fermi wavelength Af, but 
short compared to the coherence length ^o- The SN interfaces are located at 
X = zLL/2. We assume that the only scattering in the superconductors consists 
of Andreev reflection at the SN interfaces, i.e. we consider the case that the 
disorder is contained entirely within the normal region. The spatial separation 
of Andreev and normal scattering is the key simplification which allows us to 
relate the supercurrent directly to the normal-state scattering matrix. 

We first construct a basis for the scattering matrix (s-matrix). In the 
normal lead Ni the eigcnfunctions of the BdG equation (5) can be written in 
the form 



(13) 



*tc(Ni) = ( I )(fc;i)-i/2cO„(y,z)exp[±ifc°(x + iL)], 
<h(Ni)= ( ; )(fcJl)-^/^<I'„(y,^)exp[±ifc^(x+iL)], 



where the wavenumbers fcj^ and fcjj are given by 

fc^^'i^ = {2mir^f'\E^ -E^ + a-'^'sy/', (14) 

and we have defined a'^ = I, a^ = —1. The labels e and h indicate the electron 
or hole character of the wavefunction. The index n labels the modes, $n(j/, z) 
is the transverse wavefunction of the n-th mode, and En its threshold energy: 

[{pl+pl)/2m + V{y, z)]$„(y, z) - £;„$„(y, z). (15) 

The eigenfunction $„ is normalized to unity, Jdy Jdz\^n\^ = 1- The eigen- 
functions in lead N2 are chosen similarly, but with x + ■^L replaced by a; — ^L. 
In the superconducting lead Si, where A = Aq cxp(i(/)/2), the eigcnfunc- 
tions are 

V2g^)-'/'(£VAg - l)-^/^<i>„cxp[±ig,^,(x + IL)], 
1 (2gh)-V2(,2/^2 _ i)-i/4<j,^^exp[±igh(a; + U)]. (16) 



*tc(Sl) = 




ir,<=/2 
-uf/2 


<h(Si) = 


/ 
e 


i^h/2 
An^/2 


We have defined 






In"" 


^ (^ 




^c,h 


— 1 

— 2 



(2m/n2)i/2[ijp _En + n^^\e^ ^lf'^Y'\ (17) 

cr°''^arccos(£/Ao). (18) 

The square roots are to be taken such that Reg''''^ > 0, Img° > 0, Imq*^ < 0. 
The function arccosi € (0, 7r/2) for < i < 1, while arccosi = — iln[t + {t^ — 
1)^/^] for t > 1. The eigenfunctions in lead S2, where A = Aoexp(— 10/2), are 
obtained by replacing (j) by —0 and L by —L. 

The wavefunctions (13) and (16) have been normalized to carry the same 
amount of quasiparticlc current, because we want to use them as the basis for 
a unitary s-matrix. Unitarity of the s-matrix is a consequence of the fact that 
the BdG equation conserves the quasiparticlc current density 

The quasiparticlc current /qp through a lead is given by 



(20) 



One easily verifies that, provided the wavenumbers k^^ and g°'^ are real (i.e. 
for propagating modes), each of the wavefunctions (13,16) carries the same 



amount of quasiparticle current |/qp|. (The direction of the current is opposite 
for ^P+jV&j^ and ^E'^Tj^jJ".) We note that the wavefunctions (13) and (16) carry 
different amounts of charge current (i.e. electrical current), which is given by 

-^^chargo = -e— Re dy dz iu* —-u + v*—-v] . (21) 

m J J \ lox lox J 

This does not contradict the unitarity of the s-matrix constructed from these 
wavefunctions, because the BdG equation does not conserve the charge of the 
quasiparticles. (Andreev reflection [20], i.e. the reflection of an electron into 
a hole at an SN interface, is an example of a non-charge-conserving scattering 
process.) 

A wave incident on the disordered normal region is described in the basis 
(13) by a vector of coefflcients 

4" ^ {ct{N,),c-{N^),c-{N^),c+{N^)) , (22) 

as shown schematically in Fig. 1. (The mode-index n has been suppressed 
for simplicity of notation.) The reflected and transmitted wave has vector of 
coefficients 

cr ^ (c-(^i), c+(iV2), c+(iVi), c,;(7V2)) . (23) 

The s-matrix sn of the normal region relates these two vectors, c^"* — sncJJ. 
Because the normal region does not couple electrons and holes, this matrix has 
the block-diagonal form 

- ( sa{e) \ .^=( rix t 



Here so is the unitary and symmetric s-matrix associated with the single- 
electron Hamiltonian TLq. The reflection and transmission matrices r{e) and 
t{e) are N x N matrices, N{e) being the number of propagating modes at en- 
ergy e. (We assume for simplicity that the number of modes in leads A'^i and 
N2 is the same.) The dimension of SN(e) is 2A^(e) + 2N{—e). 

We will make use of two more s-matrices. For energies < £ < Aq there 
are no propagating modes in the superconducting leads Si and S2. We can 
then define an s-matrix sa for Andreev reflection at the SN interfaces by c™ = 
sac^"*- The elements of sa can be obtained by matching the wavefunctions 
(13) at |a;| = L/2 to the decaying wavefunctions (16). Since Ag ^ E-p one 
may ignore normal reflections at the SN interface and neglect the difference 
between N(s) and N(—e). This is known as the Andreev approximation [20]. 
The result is 

/ TA A ^ f e'"^/21 \ ,„,^ 



where a = exp[— iarccos(e/Ao)]. The matrices 1 and are the unit and null 
matrices, respectively. Andreev reflection transforms an electron mode into a 
hole mode, without change of mode index. The transformation is accompanied 
by a phase shift, which consists of two parts: 1. A phase shift — arccos(£/Ao) 
due to the penetration of the wavefunction into the superconductor; 2. A phase 
shift equal to plus or minus the phase of the pair potential in the superconductor 
{plus for reflection from hole to electron, minus for the reverse process). 

For e > Ao we can define the s-matrix ssns of the whole junction, by 

SSNSCg 



c2"t = ssNsc!?. The vectors 



4" ^ (c+(5i),c-(52),c;^(5i),c+(52)), (26) 

cr ^ (c-(5i),c+(52),c+(5i),c,;(52)). (27) 

are the coefficients in the expansion of the incoming and outgoing wave in 
leads Si and S2 in terms of the wavefunctions (16) (cf. Fig. 1). By matching 
the wavefunctions (13) and (16) at \x\ = L/2, we arrive after some algebra 
(using again Aq ^ E-p) at the matrix-product expression 

SSNS = U-' (1 - M)-i (1 - Mt)sNC/, (28) 

U=[ i ^^ j ,M^as^[ ^^ ^ 

One can verify that the three s-matrices defined above (sn; sa for < e < 
Aq; Ssns for e > Aq) are unitary {s^s = ss^ = 1) and satisfy the symmetry 
relation s{s, (j))ij — s(e, —(f)) ji, as required by quasiparticle-current conservation 
and by time-reversal invariance, respectively. 

We are now ready to relate the excitation spectrum of the Josephson junc- 
tion to the s-matrix of the normal region. First the discrete spectrum. The 
condition Qn == SASNQn for a bound state implies Det (1 — sasn) = 0. Using 
Eqs. (24), (25), and the folding-identity 



Detf ^ '^ j =Dct {ad- aca-'b) (29) 

(which holds for arbitrary square matrices a, 6, c, d of equal dimension, with 
Det a 7^ 0) , we find the equation 

Det [1 - a{ep)\ls„{ep)rASo{-epr] = 0, (30) 

which determines the discrete spectrum. The density of states of the continuous 
spectrum is related to ssns by the general relation [21] 

1 d 

p = —-^ TT- In Det Ssns + constant, (31) 

27ri oe 



10 



where "constant" indicates a (/)-independent term. From Eqs. (28) and (29) we 
find 

7^^ — THTTl^ilnDet [1 - a{efrlso{e)rASo{^er] , (32) 

0(p TT ocpoe 

which determines the (/)-dependence of the continuous spectrum. 

4 Short-Junction Limit 

In the short-junction hmit _L <C Ci the determinantal equations (30) and (32) 
can be simphfied further. As mentioned in the Introduction, the condition 
i <C ^ is equivalent to Aq ^ E^ , where the correlation energy Ef. = h/r^^en is 
defined in terms of the dwell time TdwcU in the junction. The elements of so(e) 
change significantly if e is changed by at least E^ [22, 23]. We are concerned 
with e of order Aq or smaller (since p{e, (f>) becomes independent of (j) for 
e 3> Ao). For Aq ^ Ec we may thus approximate so(e) ~ so(— e) ~ so(0) = sq. 
Eq. (30) may now be simplified by multiplying both sides by Det sq and using 
SqSq = 1 (unitarity plus symmetry of sq), as well as the folding identity (29). 
The result can be written in the form 



Det 



(l-£2/A2)l-ti242sin'(<^/2) =0. (33) 



For e > Aq one can see that Det[l — a(e)^r^sorASQ] is a real number. (Use that 
a is real for e > Aq and that the determinant Det (1 — aa*) is real for arbitrary 
matrix a.) Eq. (32) then reduces to dp/d(j) = 0, from which we conclude that 
the continuous spectrum does not contribute to I{(j)) in the short-junction limit. 
Eq. (33) can be solved for Ep in terms of the eigenvalues Tp {p = 1,2,... A'') 
of the hermitian N x N matrix ^12*12, 

ep^Ao[l-TpSm\4,/2)]^^\ (34) 

Since ^12^12 = ''ii4i^2i?'n {^^ follows from unitarity of so), the matrices ^12^2 
and ^21^21 have the same set of eigenvalues. We can therefore omit the indices 
of tt^ . Unitarity of so also implies that < Tp < 1 for all p. 

Substitution of Eq. (34) into Eq. (12) yields the Josephson current 

Eq. (35) holds for an arbitrary transmission matrix tt'^ , i.e. for arbitrary disorder 
potential. It is the multi-channel generalization of a formula first obtained by 



11 



Haberkorn et al. [24] (and subsequently rederived by several authors [25]-[27]) 
for the single- channel case (appropriate for a geometry such as a planar tunnel 
barrier, where the different scattering channels are uncoupled). A formula of 
similar generality for the conductance is the multi-channel Landauer formula 

[23, 28] 



„2 o„2 N 



G^^Trttt^^J^T,. (36) 

In contrast to the Landauer formula, Eq. (35) is a non-linear function of the 
transmission eigenvalues Tp. It follows that knowledge of the conductance (i.e. 
of the sum of the eigenvalues) is not sufficient to determine the supercurrent. 

5 Universal Josephson Effects 

5.1 Quantum Point Contact 

Consider the case that the weak link consists of a ballistic constriction {I ^ 
L) with a conductance quantized at G = 2NQe'^ /h (a quantum point contact 
[7]). The integer Nq is the number of occupied one-dimensional subbands (per 
spin direction) in the constriction, or alternatively the number of transverse 
modes at the Fermi level which can propagate through the constriction. Note 
that Nq <^ N. a quantum point contact is characterized by a special set of 
transmission eigenvalues, which are equal to either or 1 : 

^p ^ 1 n if Afn ^ «< Af_ ^'^'> 



ii No<p<N, 



where the eigenvalues have been ordered from large to small. We emphasize 
that Eq. (37) is valid whether the transport through the constriction is adiabatic 
or not. In the special case of adiabatic transport, the transmission eigenvalue 
Tp is equal to the transmission probability Tp of the p-th subband. In the 
absence of adiabaticity there is no simple relation between Tp and Tp. 

The discrete spectrum (34) in the short-junction limit contains an A'o-fold 
degenerate state at energy e = Aq] cos((/>/2)|. Eq. (35) for the supercurrent 
becomes 

/(0) = iVo^ sin(0/2) tanh (-^ cos(0/2)') , \cj)\ < n. (38) 

At r = the current-phase relationship is given by I{<j>) — No{eAo/h) sin((/)/2) 
for \(f>\ < TT, and continued periodically for |0| > tt. At (f> ~ tt the function 
/((/)) has a discontinuity which is smeared at finite temperatures. The ratio 
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Figure 2: Current -phase relationship at T = of a baUistic point contact 
[soUd curve, Eq. (38)], a disordered point contact (ensemble averaged) [dashed 
curve, Eq. (49)], and a tunnel junction [dotted curve, Eq. (50)]. The current is 
normalized by the normal-state conductance of the junction. 

I {(f>) / {tt AqG / e) at T = is plotted in Fig. 2 (solid curve). The critical current 
Jc = NoeAo/h is discretized in units of eAo/h. In the classical limit Nq -^ oo 
we recover the results of Kulik and Omel'yanchuk [29] for a classical ballistic 
point contact. 

Eq. (38) was derived in Rcf. [1] under the assumption of adiabatic trans- 
port. The present derivation (taken from Ref. [4]) does not assume adiabatic- 
ity. As discussed in Rcf. [1], Eq. (38) breaks down if the Fermi level lies 
within Aq of the threshold energy _E„ of a subband. In the present con- 
text this breakdown can be understood from the fact that the dwell time 
TdwoU = (i/vF)(l — En/Ep)^^^^ becomes infinitely long as Ep approaches 
En, so that the "short-junction" condition Aq ^ Ec = h/rdvjeii is violated. 
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5.2 Quantum Dot 

Consider a small confined region (of dimensions comparable to the Fermi wave- 
length), which is weakly coupled by tunnel barriers to two electron reservoirs. 
We assume that transport through this quantum dot occurs via resonant tun- 
neling through a single bound state. Let e^cs be the energy of the resonant 
level, relative to the Fermi energy Ep in the reservoirs, and let Ti/h and T2/h 
be the tunnel rates through the left and right barriers. We denote F = Fi +r2. 
If F <C AE (with AE the level spacing in the quantum dot) and T <C F/fce, 
the conductance G in the case of non-interacting electrons has the form 

where Tbw is the Breit-Wigner transmission probability at the Fermi level. 
The normal-state scattering matrix so(e) which yields this conductance has 
matrix elements [301 



lVri„F2m \ J(^Si„+S2m.) 



where X]n=i ^in = Fi, J2n=i^2n = r2. The phases (5i„, (52m, as weh as the 
basis with respect to which the matrix elements (40) are calculated, need not 
be further specified. Biittiker [31] has shown how the conductance (39) follows, 
via the Landauer formula (36), from the Breit-Wigner scattering matrix (40). 
In Ref. [3] Van Houten and the author have calculated the supercurrent 
through the quantum dot from the Breit-Wigner formula (40) in the special 
case of single- channel leads {N = 1). Just as for the conductance [31]-[33] one 
would expect that the results are not changed for multi-channel leads {N > 1). 
We will now demonstrate this explicitly in the limit F ^ Aq (which for resonant 
tunneling corresponds to the "short-junction" limit, cf. Sec. 1). The transmis- 
sion matrix product ti2t\2 (evaluated at the Fermi level e = 0) following from 
the scattering matrix (40) is 

(^124) =rBw7n7™, 7« = e'*"ViWrT- (41) 

\ / nm 

Its eigenvalues are 

rp _ j Tbw ifp=l, i.^^ 

P " \ if 2 < p < iV. ^ > 

The single non-zero eigenvalue has eigenvector Un = 7n, the other eigenvectors 
span the plane orthogonal to the vector 7. Substitution into Eq. (35) yields 
the current -phase relationship for a wide resonance (F 3> Aq), 

^ £Ao TBwsm^ ^ ^^^^^ (^^<L[i _ Tbw sin2(0/2)]V2^ . (43) 

^^' 2% [l-TBwsin2(,/,/2)]i/2 \2kBT^ ^^ ^^' >^ ) "^ > 
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Figure 3: Normalized critical current versus energy of the resonant level, at 
zero temperature and for equal tunnel barriers (Fi = r2 = Fq). The two solid 
curves are the results (44) and (45) for the two regimes Fq ^ Ag (curve a) and 
FojEros <C Ao (curve b). The dotted curve is the Breit-Wigner transmission 
probability (39). The inset shows schematically the quantum-dot Josephson 
junction. Taken from Ref. [3]. 



The critical current at zero temperature is 

eAo, 



/c = — ^[1 - (1 - Tbw)'/% if F » Ao, 
n 



(44) 



in agreement with Ref. [3]. Since Tbw = 1 on resonance (sr = 0) in the 
case of equal tunnel rates (Fi = F2), we obtain the result (3) discussed in the 
Introduction. 

For completeness, we also quote the formula for the current-phase rela- 
tionship in the opposite regime of a narrow resonance (derived in Ref. [3]): 



/c = J(4: 



ip2)i/2 [1 _ (1 _ Tbw)'/'], if F,e,es « Aq 



(45) 



As shown in Fig. 3, the lincshapcs (44) and (45) of a resonance in the 
critical current (solid curves) differ substantially from the lorentzian lineshape 
(39) of a conductance resonance (dotted curve). For Fi = F2, /c has a cusp at 
ffres — (which is rounded at finite temperatures). On resonance, the maximum 
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critical current I^'^^ equals (2eAo/?ir) min (Fi, r2) and (e/K) niin (Fi, F2) for a 
wide and narrow resonance, respectively. An analytical formula, Eq. (2), for 
the crossover between these two regimes can be obtained for the case of equal 
tunnel rates [3]. Off-resonance, /c has the lorentzian decay oc l/£ros i^i f^^ 
case F ^ Aq of a wide resonance, but a slower decay oc 1/eros in the case 
F, £res *^ Aq. Near Eros — Aq this linear decay of the narrow resonance crosses 
over to a quadratic decay (not shown in Fig. 3). 

Since we have assumed non-interacting quasiparticles, the above results 
apply to a quantum dot with a small charging energy U for double occupancy 
of the resonant state. Glazman and Matveev have studied the influence of 
Coulomb repulsion on the resonant supercurrent [2]. The influence is most 
pronounced in the case of a narrow resonance, when the critical current is sup- 
pressed by a factor F/Aq (for C/, Aq ^ F). In the case of a wide resonance, 
the Coulomb repulsion does not suppress the supercurrent, but slightly broad- 
ens the resonance by a factor ln(F/Ao) (for [/, F ^ Aq). The broadening is 
a consequence of the Kondo effect, and occurs only for Cros < 0, so that the 
resonance peak becomes somewhat asymmetric [2]. 

5.3 Disordered Point Contact 

We now turn to the regime of diffusive transport through a disordered point 
contact. We first consider the average supercurrent (averaged over an ensemble 
of impurity configurations) and then the fluctuations from the average. This 
section corrects Ref. [4]. 

5.3.1 Average Supercurrent 

The transmission eigenvalue Tp is related to a channel-dependent localization 
length (p by 

Tp^cosh-\L/Cp). (46) 

The inverse localization length is uniformly distributed between and 1/Cmin — 
l/l ior I <ti L <^ Nl [34]. One can therefore write 



dxficosh^^x), (47) 



E^=i fiTp)) /p^/^"""dx/(cosh-2 x) 



N rT.\ fi/Cmi„ , „_, -2, 



Ep=iTp) J,''-dxcosh-'x 







where (...) indicates the ensemble average and f(T) is an arbitrary function 
of the transmission eigenvalue such that f(T) -^ for T ^ 0. In the second 
equality in Eq. (47) we have used that L/Cmin — L/l 3> 1 to replace the upper 
integration limit by 00. Note that, in view of the Landaucr formula (36), the 

denominator (X]n=i -^p) i'' J^''^ h/2e'^ times the average conductance (G). 
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Combining Eqs. (47) and (35) we find for the average supercurrent the 
expression (0 < (j> < tt) 

(/(^)) = £^^(G)/°°rfx (cosh.)-sin0 



2h 2e2' 'io [l-(cosha;)-2sin2(,/,/2)]i/2 

X tanh I ^-^[1 - (cosh x) "2 sin^ (0/2)1 ^/^ 
\2kBT y 



e Jcos{4,/2) [z^ - cosmic/) /2)]y^ 

At T = this integral can be evaluated in closed form, with the result 

(lU)) = !!— ^ (G) cos(</./2) arctanh [sin(0/2)l , (49) 

e 

plotted in Fig. 2 (dashed curve). Note that the derivative dl/d(f> diverges at 
(j) = TT. The critical current is /c = max (/((/))) = 1.32 (7rAo/2e)(G), reached 
at (/) = (/)i°^ = 1.97. 

Eq. (49) for the average supercurrent in a disordered point contact agrees 
with the result of Kulik and Omel'yanchuk [35]. These authors pointed out 
that a disordered weak link and a tunnel junction with the same conductance 
exhibit nevertheless a quite different Josephson effect. For a tunnel junction 
with Tp <C 1 for all p one may linearize Eq. (35) in Tp, with the result (at zero 
temperature) 

J = ^sin<^f]Tp = ^Gsin<^, (50) 

p=i 

which is also plotted for comparison in Fig. 2 (dotted curve). The corresponding 
critical current is (7rAo/2e)G. The well-known result (50) for a tunnel junction 
[36] differs from the result (49) for a disordered weak link [35]. Zaitsev [25] has 
explained the difference in terms of different boundary conditions for the Greens 
functions at the interface between the junction and the bulk superconductor. 
The present work yields an alternative perspective on this issue: Although the 
tunnel junction and the disordered weak link may have the same conductance 
(i.e. the same value of Trtt^), the supercurrents differ because the distribution 
of transmission eigenvalues is different. In a tunnel junction all N eigenvalues 
are ^ 1, while in the disordered system a fraction l/L of the eigenvalues is of 
order unity, the remainder being exponentially small [34]. 

5.3.2 Supercurrent Fluctuations 

The analysis of Kulik and Omel'yanchuk is based on a diffusion equation for 
the ensemble- averaged Greens function, and can not therefore describe the 
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mesoscopic fluctuations of I{(t>) from the average. In contrast, Eq. (35) holds 
for a specific member of the ensemble of impurity configurations. The statistical 
properties of the transmission eigenvalues T„ in this ensemble are known from 
the theory of conductance fluctuations [11]. A general result [34, 37] is that 
a linear statistic C{f) = X]n=i fi^^n) on the eigenvalues (with f{T) -^ for 
T -^ and /(T) = 0(1) for T < 1) fluctuates with root- mean-square value of 
order unity 

rmsC={{C^)-{C)^f^^ =0{l), (51) 

independent of A^, L, or Z as long a,s I -^ L <ti Nl. An additional requirement 
is that the junction is not much wider than long: Wy, Wz ^ L where Wy and 
Wz are the widths in the y- and 2;-direction. If the junction is much wider than 
long in one direction, then rms£ is of order {WyWz/ L'^Y''^ ■ It is not necessary 
for these results that /(T) is itself a linear function of T. In the special case 
that f{T) = T, the coefficient on the right-hand-side of Eq. (51) is known [10]: 

N 

rms ^ r„ = CucF = 0.365, if Wy, Wz < L. (52) 

n=l 

According to Eq. (35) the supercurrent /(0) is a linear statistic on T„, so 
that we conclude that (at zero temperature) rms/(^) ~ eAo/^ times a function 
of (j) which is independent of the length of the junction or the degree of disorder. 
At small (/) we may linearize in T and use Eq. (52), with the result 

rms/((/)) = iCucF(eAo/?i)</'tanh(Ao/2fcBT) + 0{(ty'). (53) 

The critical current /c = max /((/)) = I{4>c) is not by definition a linear 
statistic on r„, since the phase (j)c at which the maximum supercurrent is 
reached depends itself on all the transmission eigenvalues. It is therefore not 
possible in general to write /c in the form X]n/(-^n)' ^^ required for a linear 
statistic. However, /c does become a linear statistic in the limit L/Nl <C 1 (or 
equivalently (G) 3> e^/h) of a good conductor (which is the appropriate limit 
for mesoscopic fiuctuations). To see this, we write 

I{4>) = (/(</>)) [1 + ^X{<f>)l e = L/Nl, (54) 

where the function X((j)) accounts for the sample-to-sample fluctuations of 
1(4') around the ensemble average (/(0)). One has {X) — 0, rms A — 0(1). 
We now expand /c to lowest order in e. Defining i^c = '/'c + e^c , where 
max (/((/))) = (/((/)c )), one can write 

/e = (/(0c))[l + eA(^c)] 

= (/(</.(°)+e0W)^[i + ,x(0(°)+.0W)^ 

;/(0(°))) [l + 6A(0(°)) + Oie^)] = /(0(°)). (55) 
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In the third equaUty wc have used that, by definition, d{I)/d(j) ~ at 4> — 
(j>c ■ Since I{(f>c ) is a linear statistic on r„ (note that (j)c is by definition 
independent of T„), we conclude that the critical current /c is a linear statistic 
on the transmission eigenvalues in the limit e = L/Nl — > 0, and hence that 

rms Jc = rms/((/)c ) — eAo/h. The value of the numerical coefficient remains 
to be calculated. 

Experimentally, sample-to-sample fluctuations are not easily studied. In- 
stead, fluctuations in a given sample as a function of Fermi energy Ep are 
more accessible. Josephson junctions consisting of a two-dimensional electron 
gas (2DEG) with superconducting contacts allow for variation of Ep in the 
2DEG by means of a gate voltage [38]. Point-contact junctions can be defined 
in the 2DEG either lithographically or electrostatically (using split gates) [7]. 
For such a system one would expect that if Ep is varied on the scale of Ec, the 
low-temperature critical current will fluctuate by an amount of order eAo/fi., 
independent of the properties of the junction. 
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